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Abstract 

The spin-1 Ising model with bilinear and biquadratic exchange interac- 
tions and single-ion crystal field is solved on the Bethe lattice using exact 
recursion equations. The general procedure of critical properties investigation 
is discussed and full set of phase diagrams are constructed for both positive 
and negative biquadratic couplings. A comparison with the results of other 
approximation schemes is done. 
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1 Introduction. 



The spin-1 Ising model with most general up-down symmetry, known also as a 
Blume- Emery- Griffiths (BEG) model has attracted recently a great attention as a 
simple model with rich and interesting phase structure. It was originally introduced 
Jj]] in order to explain the phase separation and superfluidity in the 3 He- 4 He mix- 
tures and latter has been developed to describe another multi-component physical 
systems, such as methamagnets, liquid crystal mixtures, microemulsions, semicon- 
ductor allows, etc. 

The model is defined by Hamiltonian: 

-m = JJ2 SiSj+KYl ^ 2 -a$> 4 2 (i) 

<ij> <ij> i 

where takes the values ±1, at each lattice site, (ij) denotes a summation over all 
nearest-neighbour pairs, J and K are correspondingly the bilinear and biquadratic 
interaction constants, A is a single-ion crystal field. 

The spin-1 Ising model was solved exactly only on two-dimensional honeycomb 
lattice in a subspace of interacting constants e K cosh J = 1 [^, |3], |j] (recently such 
solution has been found for higher spin-S models as well [§, |], [7|]) but its critical 
properties for positive J, K > were well established by different approximation 
techniques Jl|, || ^ [10| . The most interesting result here was the first occurrence (in 
theoretical model) of the tricritical point, at which the second order phase transition 
line (A-line) turns to the first order one. 

The picture is practically not changed for negative bilinear couplings J < 0. On 
the bipartite lattice (i.e. the lattice which can be divided on two sublattices A and 
B, such that every site belonging to A is surrounded only by sites belonging to B 
and vice versa) the region J < is mapped on the region J > by redefining the 
spin directions on one sublattice. As a result we obtain the same phase diagrams, 
where only ferromagnetic phase is replaced by antiferromagnetic one. Hence with 
no loss of generality we can consider only the case J > 0. 

On the other hand, the negative values of biquadratic coupling K change the 
situation drastically. The region Kj J < is now a subject of intensive studying. A 
new staggered quadrupollar phase (also called antiquadrupolar) was predicted and 
investigated on the square lattice by means of mean-field approximation (MFA) and 
by Monte Carlo (MC) simulations The direct first order transition was found 



from antiquadrupolar (a) to ferromagnetic (f) phase, but it was not confirmed by 
the recent MC (jl2[ and cluster variation method (CVM) |y| studies. Now it seems 
that a and / phases are always separated in two dimension by disordered phase (d) 
and they meet only at T = 0. This direct a <-> / transition was however established 
on three-dimensional cubic lattice [[TJ], [T5], M. The global MFA analysis [TJJ] on this 
lattice showed also a number of other remarkable features, such as doubly reentrant 
behaviour at > Kj J > — 1 and new staggered ferrimagnetic phase which appears 
between a and / phases at K/J < —1. Latter investigations J0| |17|, [18], 0, ^(J 
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mainly confirmed these results, however a number of contradictions still remain. 
This makes interesting the further consideration of the model, especially by applying 
another approximation tools. 

In the present paper we consider the solution of the spin-1 Ising model on the 
Bethe lattice. We review both positive and negative values of biquadratic coupling 
K and show that, though in main details we reproduce the phase diagrams obtained 
by other authors, there are some essential differences, concerning the place and order 
of phase transitions. 

The paper is organized as follows. In the Sect. 2 we introduce the model and de- 
rive some analytical expressions including the set of exact recursion equations. The 
procedure of critical properties investigation based on these equations is described 
in Sect. 3. Resulting phase diagrams are presented and discussed in Sect. 4. Final 
Sect. 5. is devoted to Conclusion. 



2 Model formulation. 

The BEG model is characterized by two order parameters, magnetization m and 
qudrupolar moment q: 

m = («*>, q = (a?) (2) 

However to account the possible two-sublattice structure we need actually four 
order parameters: m^ B = (sj)a,b and q^B = (s^)a,b, where A, B denotes sublat- 
tices. These parameters define the four different phases of the BEG model: 

1. disordered phase (d) : rriA = wig = 0, qA = Qb 

2. ferromagnetic phase (/) : rriA = tub ^ 0, Qa — Qb 

3. antiquadrupolar phase (a) : rriA = ttib = 0, g^/ Qb 

4. ferrimagnetic phase (z) : ^ rriA ^ rriB ^ 0, qA ^ Qb 

The Bethe lattice consideration for any model is based, in some way, on one or 
more exact recursion equations. We construct these equations in the following way 
[PH| . Taking into account the shell structure of Bethe lattice (Fig. 1) one can express 
the partition function Z of the model on the finite n-shell lattice in the form: 

Z = £ -m = £ex P (-A S ^ n ( So )r (3) 

{s} s 

where z is a lattice coordination number, sq denotes the central spin and g n (so) is a 
contribution to the partition function of one lattice brunch, starting from the central 
site with fixed spin value s . The latter is readily connected with gr n _ 1 (si) : 

9n(s ) =5Zexp(Js si + ATsgs?- As?)[p n _i(si)] 2_1 (4) 

si 
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Introducing new notations: 



_ 9n{ + > 9n 



9n(0) 



9n(0) 



(5) 



and summing up over all values of central spin sq (i.e. ±1, 0) we obtain a set of two 
recursion equations: 



ip{x 

ni Un) 



where 



Vn+l = <P(Vn, X n ) 

e A + e K (e J u z ~ 1 + e _ V -1 ) 



(6) 



<p(u,v) 



e A + u z ~ x + v z ~ x 



The values x and y have no direct physical sense, but one can express in terms 
of x and y all thermodynamic functions of interest. Thus order parameters (0) will 
be written in the form: 



m 



x z — y z 
e A + x z + y z 

x z + y z 
e A + x z + y z 



Using equations we can writte also the expression for the free energy: 

-Pf-jfhZ 

in the form: 

-Pf = ]n [l + e- A (x^ + ^)l + 



2-z 



In 



1 + e-^x*- 1 + y'- 1 ) 



(7) 
(8) 

(9) 
(10) 



3 Critical properties investigation. 

The equations @ form an iteration sequence {x n ,y n }, which in the thermody- 
namic limit converges to stable fixed points. Via the expressions (|7-10) these 
points completely define the possible states of the system. The remarkable points of 
this approach is that non-staggered phases are described by the single fixed points 
{x n ,y n } — ► {x,y}, while the staggered phases appear as a 2-cycle doubled points 



{%m yn} 



{x A ,y A } for odd n 
{xb, for even n 



This property can be explained by the fact that all sites of each individual shell 
of the Bethe lattice belong to the same sublattice. 
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Thus, noting from (|7|) that m = means x = y, we can, in our case, classify the 
four mentioned phases as follows: 



1. disordered phase (d) : x = y, single fixed point; 

2. ferromagnetic phase (/) : x ^ y, single fixed point; 

3. antiquadrupolar phase (a) : x — y, period doubling; 

4. ferrimagnetic phase (i) : x ^ y, period doubling. 

It is possible to obtain the full bifurcation picture, including chaos on some 
hierarchical lattices [^|, but not in the bipartite case (which is the Bethe lattice 
with nearest-neighbour interactions). We have only first period doubling and it is 
well known from the theory of iteration processes [^4 that all our points of interest 



(i.e. stable fixed points and 2-cycle doubled points) can be found among the solutions 
of the set of four equations: 

x A = ip(x B , y B , J, K, A) (i) 

Va = <p(vb, %b, J, K, A) (ii) 

x B = ip(x A , Va, J, K, A) (iii) 

Vb = viVA, xa, J, K, A) (iv) 



The physical stable solutions of this set define the pure states of the model. As 
a matter of fact there is no need to solve ([U]) in general form. Knowing apriori the 
possible phases of the model we can separate the solutions of (|TT| ) concerning the 
given phase. Thus the disordered phase (x — y, A = B) can be defined by single 
equation, either (i) or (ii). This phase is not degenerate. For ferromagnetic phase 
(x 7^ y, A = B) it is enough to consider two equations, namely (i) and (ii) and 
to exclude the disordered solution. This phase is doubly-degenerate (m <-» —to) 
due to obvious symmetry of these two equations under the (x <-> y) transformation. 
The antiquadrupolar phase (x = y, A ^ B) again can be found from equations (i) 
and (iii), after excluding the disordered solution. This phase is two-fold degenerate 
because of the A <-> B symmetry. And only for ferrimagnetic phase we have to 
consider all four equations, but they are simplified by excluding previous solutions. 
Ferrimagnetic phase is four-fold degenerate since the (x <-> y) and (A <-> B) sym- 
metry. Of course, the last two phases are also infinitely degenerate with non-zero 
residual entropy. 

Further procedure can be roughly described by the following steps: 

1) The intersections of the solutions, describing the different phases, give us, 
generally speaking, the points of second order transitions. 

2) The presence of several simultaneous solutions at the given K, J, A again 
generally speaking means the co-existing phases and first order transition, which 
should be located by matching free energies flit]) of these phases. 

3) The intersections of first and second order critical lines give the critical and 
multicritical points of several types. 
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On Fig. 2 we illustrate this procedure on a simple example of appearance of the 
tricritical point. The disordered and ferromagnetic solutions of (|TT|) are shown on 
a crystal field A versus quadrupolar moment q plot for different fixed temperatures 
1/zJ. In high temperature region the d and / solutions intersect in their stable parts 
and this leads to a second order transition. At a low temperature the d line does 
not meet the stable part of / line and first order transition takes place. Tricritical 
point is observed in an intermediate situation. 



4 Phase diagrams 

The resulting phase diagrams for spin-1 Ising model are constructed on the Bethe 
lattice with coordination number z — 4. As it is common we plot them as a distinct 
constant Kj J cross-sections, in the temperature 1/zJ versus crystal field A/zJ. We 
obtain eight qualitatively different diagrams in the whole range of parameter space 
(Fig. 3). 

The > K/J > — 1 counterpart is characterized by absence of doubled points of 
the recursion sequence {x n , y n } and hence represents the non-staggered region of the 
BEG model with two phases d and /. The first three phase diagrams (Fig. 3a-c) are 
well known from the pioneer work of Blume, Emery and Griffiths |J. Our results 
are quite agree with general picture. The tricritical point doesn't appear for large 
positive K/J (Fig. 3a), since the second order line, limiting ferromagnetic phase 
from above, terminates earlier at the critical point E by the first order line limiting 
/ phase from the right. The latter line itself terminates at critical point C, and in 
higher-temperature segments two subphases of disordered phase with different dense 
co-exist. For the K/J close to 3 (Fig. 3b) we observe tricritical point T at which 
second order line (A-line) turns to the first order one and also the triple point R, 
where three different first order transitions meet. This structure disappears as K/J 
goes down, and simple tricritical point is seen (Fig. 3c) for Kj J close to (from 
both positive and negative sides). 

Very interesting critical phenomenon, called doubly reentrant behaviour, takes 
place for the values —0.35 > Kj J > — 1 (Fig. 3d). At a fixed crystal field AjzJ the 
model exhibits the disorder-ferromagnetic-disorder-ferromagnetic sequence of phases 
as temperature is lowered. The dependence of order parameters versus temperature 
is shown for this region on Fig. 4. Doubly reentrant structure shrinks to a zero 
temperature and at K/J = — 1 the only second order line remains, which reaches 
the T = axes at a point A = (Fig. 3e). 

Comparing our results with other approximations, we note the following differ- 
ences. Doubly reentrance appears in our model for lower values of Kj J and continue 
till Kj J = — 1, hence we do not observe the internal critical point structure of fer- 
romagnetic phase. This structure was found by MFA fTjl and confirmed by RG 



studies [fL8| , but it was also not established by CVM [20]. Besides we would like to 
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mention that nowhere in region > K/J > — 1 we obtain the single reentrant part, 
though it was found by other authors. 

The 2-cycle doubled solutions of flTTD appear at K/J < —1, and this means the 
emergence of new phases with broken sublattice symmetry. 

At —3 < K/J < — 1 (Fig. 3f) the antiquadrupolar phase is separated from d 
phase by second order line and from / phase by first order one. These two lines 
meet with the second order line, separating the d and / phases, at bicritical point 
B. The ferrimagnetic phase lies at the low-temperature region and is separated by 
the second order line from a phase and by first order line from / phase. These two 
lines meet with a-f first order line at critical end point E and at point S. The 
latter point locates at T = 0, A/zJ = K/J + 1 and describes the macroscopically 
degenerated ground state with non-zero residual entropy. 

As Kj J lowers the point E approaches to the B point and they coincide at 
K/ J = —3 in the point A (Fig. 3g), such that there is no direct transition from a to 
/ phase. Thus A is a new multicritical point at which three second order lines and 
one first order line meet. Note that the first order i-f transitions locate at straight 
vertical line and this locus K/J = —3, A/qJ = —2 corresponds to zero-field 3-state 
antiferromagnetic Potts model with Hamiltonian: 

-/3W = -2J £ 5 Si7Sj , J>0 (12) 

<ij> 

At K/J < —3 the transitions from i to / phase in high temperature region 
become of the second order (Fig. 3h). As a result the new tricritical point T' arises 
inside the ordered region. The multicritical point A turns to the tetracritical point 
M, at which four second order lines meet with different slopes. 

Though in general the last three diagrams are similar to those obtained by MFA, 
there is one essential difference: in our approach the transitions from a to i phase are 
of second order (first order in MFA), while the transitions from i to / phase are of first 
order at least for Kj J > —3 (always second order in MFA). This means in particular 
that we observe another type of ferrimagnetic phase, which co-exists with / phase 
and is caused by instability of a phase against the spontaneous magnetization. Such 
phase was found by CVM |0], P0f| , but only in high-temperature region, while it is 
the only ferrimagnetic phase present in our consideration. We would like also to 
mention the very narrow region of occurrence of ferrimagnetic phase. 



5 Conclusion 

Using exact solution on the Bethe lattice we have constructed the full set of phase 
diagrams for spin-1 Ising model for both positive and negative biquadratic coupling 
K. These diagrams feature all recently found properties of the model, including 
doubly-reentrant behaviour, staggered quadrupolar and ferrimagnetic phases and 
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great number of different critical and multicritical points. Thus we can say about 
the validity of Bethe lattice approximation and summaries some advantages of this 
method. The quantitative comparison of the results has not been our purpose, but 
Bethe lattice solution was shown to be more exact than MFA f25| , p6| . Besides it 
is quite easy to use and it provides the analytical expression of all thermodynamic 
functions of interest, so the complete information about the system under the study 
can be obtained. 

As to the mentioned disagreements, especially in the staggered region, they may 
be caused by dimensionality effects JL5|]. The Bethe lattice is effectively infinite 
dimensional, but one can successfully approximate the real lattices in different di- 
mensions by changing the coordination number (this was shown in particular in 
the global consideration of the antiferromagnetic Potts model P7fl). We have con- 
structed the phase diagrams for the Bethe lattice with coordination number z = 4. 
We do not consider the simplest case z = 3 since it leads to qualitatively different 
phase diagrams, similar to those, which were obtained on two-dimensional lattices 
(see Introduction). Our preliminary study shows also that phase diagrams change 
in case of greater coordination number. This question we are going to clarify in our 
further work. 
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Figures captions: 



Fig. 1. Bethe lattice with coordination number z = 3 

Fig. 2. Graphical representation of the solutions of equation (|TT|) on a A versus q 
plot at different temperatures: 1/zJ = 0.4 (a), 1/zJ = 0.15 (b), 1/zJ — 0.23 
(c) and constant K/J — 0. d and / are respectively the disordered and 
ferromagnetic phases. Solid lines represent the stable solutions, dashed lines - 
unstable and non-physical solutions. Thick solid lines show the main behaviour 
of the system. 

Fig. 3. Phase diagrams of the spin-1 Ising model on the Bethe lattice with coordi- 
nation number z = 4 at constant K/J values: 5 (a), 3 (b), -0.1 (c), -0.8 (d), -1 
(e), -2.5 (f), -3 (g), -3.5 (h). Disordered d, ferromagnetic f, antiquadrupolar a 
and ferrimagnetic i phases are present. Dashed and solid lines indicate respec- 
tively the first-, second order transitions and C, E, R, T, T'S, B, A, M indicate 
the critical and multicritical points of different types (see in text). Some fair 
details are shown in the insets. 

Fig. 4. The dependence of order parameters m,q versus temperature 1/zJ at a 
constant A/zJ = 0.099, K/J = —0.8 in the doubly reentrant region. Close, 
open arrows indicate respectively the places of first-, second order transitions. 
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